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possibility of controllable generation of given quantum states. In this paper, we demonstrate creation 
and detection of specific ground-state magnetic-sublevel superpositions in 87 Rb. By appropriate 
tuning of the modulation frequency and magnetic-field induction the efficiency of a given coherence 
generation is controlled. The processes are analyzed versus different experimental parameters. 
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I. INTRODUCTION 

Quantum coherences are at the heart of many nonlin- 
ear and quantum optical phenomena. They are respon- 
sible for the appearance of such effects as coherent pop- 
ulation trapping [l(, electromagnetically induced trans- 
parency (3 , extremely slow light propagation Q and its 
storage [4] in a medium, etc. A possibility of control- 
lable generation and modification of an arbitrary quan- 
tum state also lays at the very foundations of quantum- 
state engineering and quantum-information processing 
(see, for example, Ref. |5|). 

A specific example of quantum coherence phenomenon 
is nonlinear magneto-optical rotation (NMOR) @, 0|- 
The effect consists of light-intensity-dependent rotation 
of the polarization plane of linearly polarized light upon 
its propagation through a medium placed in an external 
magnetic field. The effect is based on the generation, evo- 
lution, and detection of non-equilibrium population dis- 
tribution and/or quantum coherences between Zeeman 
sublevels of a given atomic state. In a typical Faraday 
geometry, in which the magnetic field and light propaga- 
tion direction are parallel, linearly polarized light gener- 
ates coherences between Zeeman sublevels differing in the 
magnetic quantum number m by even values, Am = 2n,i, 
where rii is an integer [8|. Despite the fact that differ- 
ent types of coherences can be generated in atoms, it is 
not usually possible to separate contributions from co- 
herences with particular Am to the NMOR signal that 
is observed around zero magnetic field, B m |12l. [l3|. 

An important breakthrough in the study of NMOR 
was the application of frequency- (l5| and amplitude- 
modulated light which resulted in the FM NMOR 
(frequency modulated nonlinear magneto-optical rota- 
tion) and AMOR (amplitude modulated optical rota- 
tion) techniques. Application of the modulation tech- 
nique enables the generation of given types of atomic 
coherences, i.e., coherences between Zeeman sublevels 
with specific m. By exploiting spatial symmetries of the 
atomic angular-momentum distribution associated with 
a given quantum state, it is possible to selectively cre- 
ate superpositions between sublevels differing in a mag- 
netic quantum number m by 2, 4, or even more, if only 
the system supports such coherences |H, [l(| . Information 



about the system's quantum state may also be obtained 
from the angular momentum distribution by analyzing 
time-dependent rotation of the polarization plane of light 
measured at non-zero magnetic field. Thus, in addition 
to selective generation and detection of the coherences, 
the technique constitutes a powerful tool in analyzing the 
evolution of a quantum state of a system. In particular, it 
allowed detailed investigations of relaxation processes of 
ground-state coherences in atoms contained in a paraffin- 
coated vapor cell H2HHJ. 

This article presents the investigations on the genera- 
tion, evolution, and detection of long-living ground-state 
observables (non-equilibrium population distribution and 
ground-state coherences, the latter represented by non- 
diagonal density matrix elements between the ground- 
state sublevels) in atomic vapor subjected to a longitu- 
dinal magnetic field. The measurements are performed 
in the pump-probe arrangement with one beam used for 
the creation of atomic polarization and another beam em- 
ployed for the detection of the system's quantum state. 
With this arrangement we create the Am = 2 Zeeman 
coherences in non-zero magnetic fields. By appropriate 
tuning of the pumping laser, the coherences are gener- 
ated in the F = 1 and F = 2 hypcrhne ground states 
of 87 Rb. These coherences evolve in the external mag- 
netic field with the frequency determined by the energy 
splitting between respective sublevels (the Larmor fre- 
quency), and are continuously probed with CW light. 
We study the generation and evolution of the coherences 
at various pumping and probing conditions. Using mag- 
netic fields such that the nonlinear- Zeeman splitting of 
the ground-state sublevels is comparable to or exceeds 
the relaxation rate of this state coherences, we resolve 
and selectively analyze all three Am = 2 coherences gen- 
erated in the F = 2 state. 

The article is organized as follows. The next 
section presents a theoretical approach demonstrating 
the relation between ground-state Zeeman coherences 
and nonlinear magneto-optical rotation with amplitude- 
modulated light. Section IIIII describes experimental ap- 
paratus, while Sec. HVI presents the results and their anal- 
ysis. Final remarks are collected in Sec. [V] 
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II. THEORETICAL RELATION BETWEEN 
AMOR SIGNAL AND GROUND-STATE 
COHERENCES 

The optical properties of a medium are characterized 
by the complex refractive index rj 

T] = n + in= 1/I + x, (1) 

where n, k, x denote, respectively, the refractive index, 
absorption coefficient, and medium electric susceptibil- 
ity. Knowledge of r\ for different polarizations allows one 
to determine anisotropic properties of the medium. In 
particular, the Faraday angle tp, which is the angle of 
polarization rotation upon transition through a medium, 
may be calculated using refractive indices n + and n_ of 
two circular polarizations er + and a~ 

¥ = i^(n+-n-), (2) 
I 



where u> is the light frequency, L the length of the 
medium, and c the speed of light. In order to calculate 
the susceptibility of a medium, and hence the refractive 
indices and Faraday angle, the density-matrix formalism 
may be used 



X± — ~pr ), d mm ±ip m ±i m , (3) 

u m =-F 

where p m ±\ m is the optical coherence between the 
ground-state sublevel to and excited state sublevel m ± 1 
and d mm ±i is the corresponding dipole matrix element, 
iV the number density of atoms, and Eq is the amplitude 
of the electric field of the light. Substituting Eq. ^ into 
Eq. (TTJ) and expanding the equation into the power series 
with respect to x allows one to link the Faraday angle 
with the density matrix elements 



f = 



uLN 
2cE 



F-l 



Re d 



mm-\-l Pm+lm 



\?n — — F 



- E d 

m=-F+l 



mm— lpm— lm 



(4) 



r 



In order to calculate the density matrix elements, the 
evolution of the density matrix p, governed by the Liou- 
ville equation, needs to be considered 



Q 



- h [H, g] 



A, 



(5) 



where H is the total Hamiltonian of the system, T the 
relaxation operator, A the repopulation operator describ- 
ing p-independent mechanisms such as transit and wall 
relaxation |2Jj] , and the square and curl y b rackets denote 
the commutator and anticommutator [2l| . In the con- 
sidered case, the total Hamiltonian of the system is a 
sum of the unperturbed Hamiltonian Hq and the Hamil- 
tonians describing interactions of atoms with light Vi and 
magnetic field Vb 



H = Ho + Vi+Vb. 



(G) 



Assuming the y-polarized laser beam, the Hamiltonian 
describing the light-atom interaction may be written in 
the dipole approximation as 

Vi = -E • d = -E a e- %ult d y = — ]=E ( 



v/2 



(7) 

where E is the electric field of the light of amplitude 
Eq, d denotes the electric dipole moment operator, and 
d± are the dipole-matrix elements corresponding to the 
transitions between the ground- and excited-state Zee- 
man sublevels differing in the magnetic quantum number 
to by ±1. 



The magnetic-field interaction Hamiltonian Vb may be 
presented in the form 



V B = n ■ B, 



(8) 



where \x is the magnetic dipole moment operator, and 
B the magnetic-field induction. Since in our geometry 
the magnetic field is applied along the quantization axis, 
it changes the Zeeman-sublevel energies removing their 
degeneracy but does not mix them. For an alkali atom, 
one can calculate the energy shift b of a given ground- 
state magnetic sublevel to p using the Brcit-Rabi formula 



LO B (rn F ) = — 
n 



2(27 + 1) 



± 



4m 
27+1' 



(9) 

where x = uj^/Abf with tu^ = gppBB/Ti being the Lar- 
mor frequency, gp the Lande factor of the state with 
a total angular momentum F, ps the Bohr magneton, 
Ahf the energy splitting of the ground-state hyperfine 
doublet, E u the "center of mass" energy of the level with 
no hyperfine interaction, and the signs ± correspond to 
two hyperfine components F = I ± 1/2. 

Combining Eqs. ([H])-© with Eq. allows one to for- 
mulate equations describing time evolution of a given 
density-matrix element p a b 

Pab = -i^abPab+i ^ i^ajPjb ~ Paj^jb)-^ab (Pab ~ pH) , 

(10) 
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where Q a j — E a jd a j /y/2h is the Rabi frequency asso- 
ciated with the transition between \a) and |j) states, 
uj ab = ujB(m a ) — ws(mb) denotes the frequency splitting 
of the levels, and p^? the equilibrium value of p ab . 

In order to demonstrate the role of ground-state coher- 
ences in rotation of the polarization plane, the explicit 
formulae for optical coherences need to be written. To 
derive such analytical formulae, we use the perturbation 
approach, where the density matrix is expanded into the 
power series of the electric field of the light 

oo 

p = Y J p {n] K- (ii) 

71=0 

In such a case, the relation for time evolution of a given 
density matrix element takes the form 

Pab = -^abP ab ' + I {Majp) b - P a] fyb ) - 
3 

r ( (n) eg\ 
~ T ab [Pab - Pal) ■ 

(12) 

Even though in this paper nonlinear magneto-optical 
rotation with amplitude-modulated light is studied, we 
first derive analytical formulae for populations, optical 
and Zeeman coherences when CW light is used (NMOR) . 
Such an approach facilitates the understanding of the 
problem and provides an insight into the modulated case. 
By simple generalization of the CW formulae one can 
write relations for the density-matrix elements when AM 
light is applied. 

A. Low magnetic field, unmodulated light 

Application of the unmodulated light allows one to em- 
ploy the steady-state approximation (p = 0) . Within this 
approximation the density-matrix elements calculated in 
the first three orders of the expansion are given by 

= p eq 

raa raai 

(1) = n ab (0) 

ab A , , , ."p raa ' 

n (!) Wo 

(2) _ "a&gfcq - V ab "fca 

PaCl Taa ' (13) 

O (!) Wo 

(2) = U abV ba , - O ah Uba' 

aa ' —i i , _l_ f P , ' 

uJaa' t t-i- aa' 

O ( 2 ) i ( 2 )o 

(3) _ MabPaa + P aa '* l a'b 

° ab ~ Auj ab + tT ab ' 

where p aa denotes the population of the ground-state 
sublevel \a), a ab the amplitude of the optical coherence 
Pab (Pab = cr a be~ lult ), p aa ' the ground-state coherence, 
Au> ab = uj — uj ab the light detuning from the transition be- 
tween the ground state \a) and the excited state |6), and 
the superscript to p denotes the order of the expansion 



[22l |. From Eqs. (|13l) . the amplitude of the third-order 

(3) 

optical-coherence a ab depends on the ground-state Zee- 

man coherence p ', which is characterized by the ground- 
state relaxation rate r aa '- Since the ground-state relax- 
ation is much slower than the relaxation of the optical 
coherence, T aa i <C r ao , the third-order optical coherence 
manifests in the absorption and dispersion via spectral 
features much narrower than those associated with the 
first-order coherence. Moreover, Eqs. fTB")) additionally 
show that not only the widths but also the light-intensity 
dependences differentiate between the two contributions; 
while the first-order optical coherence is linear in J7, and 
thus is the amplitude of the electric field Eq, the third- 
order coherence depends on fi 3 . Thus, based on the 
intensity dependences, one demonstrates that the first- 
order optical coherence is responsible for linear optical 
phenomena, such as polarization rotation independent of 
light intensity, and the third-order coherence determines 
nonlinear phenomena like NMOR and EIT. 

As shown above, NMOR is associated with ground- 
state Zeeman coherences. In particular for the F — 
2 state, three different ground-state coherences (p-2,0, 
P-11, Po,2) contribute to NMOR. Thus, the effect may be 
used for investigation of the ground-state coherences, in 
particular, their generation and evolution under interac- 
tion with external fields. It is noteworthy, however, that 
at low magnetic fields independent studies of a given co- 
herence are not possible because of the same dependence 
of all contributions on the magnetic field and light inten- 
sity. Such a distinction would be possible for higher mag- 
netic fields when Zeeman sublevels depend nonlinearly on 
a magnetic field, however, at such fields the NMOR sig- 
nals are not observed. 



B. Stronger magnetic field, modulated light 

When the intensity of light is sinusoidally modulated, 
the electric field of the light takes the form 



E= — -j= — VI - a m cos uj m t, (14) 

where a m is the modulation amplitude and uj m the mod- 
ulation frequency. It may be easily shown that for the 
full modulation (a m = 1), relation (|14|) simplifies to 
E = Eocxp(—iu}t)sm(u} m t/2) This modification of 
the light spectrum results in a change of the light-atom 
interaction Hamiltonian 

Vi=- ( e - i{ " + " m/2)t + e -<»-»™/ 2 )A (d_ + d+) . (15) 

Application of the modulated light rules out the stan- 
dard steady-state approximation. This is caused by 
the appearance of the time-dependent Rabi frequency 
n = Qe lu ™/ 2t + Q e - lul ™ /2t that drives oscillation of the 
density-matrix elements at different frequencies. In or- 
der to solve Eq. (fTU)) , the density matrix needs to be 
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expanded into the Fourier series of the modulation fre- 
quency u) m /2 

OO 

P= E p (k) e lk ^' 2t , (16) 

k— — oo 

where p^ is the fc-th Fourier coefficient. Introduction 



of the Fourier expansion (JTHJ) into Eqs. (IT5)l enables ap- 
plication of the steady-state approximation for a given 
Fourier coefficient of the density matrix. In such a case, 
one can calculate the time-dependent density matrix ele- 
ments j o" ,fe ) [superscripts (I, fc) denote the l-th. order of 
the perturbation expansion and the fc-th order of the 
Fourier expansion in half the modulation frequency] 



p (0,k) 
raa 



p (2,fc) 
raa 



(2,fc) 
Paa> 



T (3,fc) 
7 ab 



Pal S k,0: 



Qal 



(0,fc-l) (0,fc+l) 
Paa T yaa 

Aw ab - kuj m /2 + iT ab ' 
ba 



n ( (l,fc-l) . (l,fc+l)\ n / (l.fc-l) . (i,fc+i)\ 
Va b [<7 ba + K a ') - Mba [Kb + Kb ) 

-ku) m /2 + iT aa 
n ( (l.fc-i) . n / (i,fc-i) (i,fc+i)\ 

^ab [K a , + K a , ') - Mba' [Kb -Kb ) 



(17) 



-UJ aa ' - kui m /2 + iT aa , 
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, (J?*-l) (2,fe+l) 
ab 1 paa i f J aa 



(!) ( J 2 ^ 1 ) , n ( 2 . fc + 1 ) 
"a'b I Paa' ' Paa' 



Au ab - ku m /2 + iT ab 



where Si m is the Kronecker delta. Although relations for 
CW and AM light [Eqs. (TT3]) and (fT7])] are similar there 
are some significant differences between them. First is 
the appearance of cross terms that couple different or- 
ders of the Fourier expansion. For instance, the fc-th 
order populations and Zeeman coherences couple to the 
fc ± 1-th orders of optical coherences. Since the largest 
density-matrix elements are those with low fc (in zeroth 
order the only non-zero elements are the ground-state 
populations), the coupling to the higher-order density- 
matrix elements is weaker. This enables truncation of 
the formally infinite series (fT6|) at some finite fc c (usually 
not bigger than 5). The fc-(fc±l) dependence additionally 
results in the zeroing of some density matrix elements. It 
may be shown that populations and Zeeman coherences 
are nonzero only at even fc (for odd fc p^ffi — 0) and the 
only non-zero optical coherences are these evaluated at 
odd fc. The last difference manifests in the appearance 
of the — fcwm/2 term in denominators of the formulae for 
the density-matrix elements, which leads to the genera- 
tion of additional resonances of the density-matrix ele- 
ments vs. the modulation frequency. Assuming that the 
whole energy splitting of the ground-state Zeeman sub- 
levels is due to the magnetic field the resonance arises 
at non-zero magnetic field; the time-dependent rotation 
of the polarization plane arises when the Larmor split- 
ting of the levels coincides with a given multiplicity of 
half of the modulation frequency (fcw m /2 = — cj Qa /). For 
F > 1 and low magnetic field all Am = 2 ground-state 
coherences are generated with equal efficiency (same en- 
ergy splitting of the levels) and single AMOR resonance 



is observed. At stronger fields, i.e., when the nonlinear 
Zeeman splitting of the sublevels is comparable to, or 
exceeds the ground-state relaxation rate, each Am = 2 
coherence has different resonance frequency. It is this 
difference, which allows selective addressing of a given 
ground-state Zeeman coherence. 

In order to calculate dynamic nonlinear magneto- 
optical rotation, i.e., the AMOR signal, one needs to 
combine all Fourier coefficients of the density matrix and 
introduce them into Eq. Q 



x±(t) = 



(18) 



E ( e -iu m /2t _|_ e iu m /2t\ 

Multiplication of Eq. (|18D by sin(mo;„ i i) or cos(mu; m i) 
and integrating them over the modulation period gives 

X±(t) sin(mw TO £), 

(19) 



(m) 
X± An 



(rn) 
quad 



2tt /mujj, 



X±(t) cos(mw m i), 



where m is the harmonic number, allows one to find for- 
mulae for the in-phase/ quadrature amplitude of electric 
susceptibility at a given harmonic of the modulation fre- 
quency. Substituting Eqs. (fT9")) into Eq. ([T]) and then 
into Eq. (01 enables calculation of the amplitude of time- 
dependent nonlinear magneto-optical rotation (AMOR 
signal). It should be noted that the AMOR signal mea- 
sured in our experiment, i.e. at the first harmonic of 
the modulation frequency (rn — 1), is described by the 
third-order optical coherence and fc = 1. 
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III. EXPERIMENTAL APPARATUS 

The layout of the experimental apparatus is shown 
in Fig. [U A paraffin-coated buffer-gas-free cylindrical 
glass cell of 2-cm in diameter and length of 2 cm con- 
tains isotopically enriched sample of 87 Rb. The cell is 
heated to 50° C by a non-magnetic resistive oven provid- 
ing atomic density of about 5 x 10 10 atoms/cm 3 [25j | . 
The cell is placed inside a three-layer mu-metal mag- 
netic shield reducing the external, uncontrollable mag- 
netic fields by a factor higher than 10 4 . The resid- 
ual fields are compensated with two sets of orthogonal 
magnetic-field coils: one for the first-order magnetic-field 
gradients and another for the second-order gradients. An 
additional solenoid is used to generate a highly homoge- 
nous and well controlled magnetic field along the light 
propagation direction which is varied within a range of 
±1 G. 

The rubidium atoms interact with two co-propagating, 
linearly polarized light beams: the pump and the probe. 
Both beams are generated with the same external-cavity 
diode laser but their intensities are controlled indepen- 
dently. The laser-light frequency is monitored with 
a saturated-absorption-spectroscopy system and can be 
stabilized to a particular transition of the Rb Dl line 
(795 nm) with a dichroic-atomic- vapor laser lock [29| . 
The intensity of the pump light is modulated with a 
single-pass acousto-optical modulator (AOM) optimized 
for the first order diffraction. Application of AOM en- 
ables modulation of light with an arbitrary frequency, 
amplitude, and waveform. It also leads to a frequency 
shift of the pump light relative to the probe by 80 MHz. 
After traversing AOM, the pump light illuminates atoms 
contained in the vapor cell. The atoms are simultane- 
ously probed with the unmodulated light beam, split off 
from the main beam before AOM. A balanced polarime- 
ter situated after the shield is employed to analyze the 
polarization state of the probe. A small angle between 
the beams allows blocking of the pump before the po- 
larimeter. The polarimeter consists of a Glan polarizer 
and two photodiodes. The polarimeter differential sig- 
nal is demodulated with a lock-in amplifier at the first 
harmonic of the modulation frequency. This signal is 
than electronically divided by twice the sum of photodi- 
ode signals which, for not-too-big rotations, yields infor- 
mation about the amplitude of the polarization rotation 
[tp rj (7i — I%)/2{Ii + I2), where Jx,2 are the respective 
light intensities in the first and second channel of the 
polarimeter]. Eventually, the signal is stored on a com- 
puter which also controls the light modulation and the 
magnetic-field strength. 



IV. RESULTS AND DISCUSSION 

Figure [5] shows the amplitude and phase of the AMOR 
signal measured vs. the modulation frequency for the 
pump tuned to the center of the F — 2 — > F' = I 




FIG. 1: Experimental setup. D is the detector, IO the optical 
isolator, AOM the acousto-optical modulator, P the polar- 
izer, BS the beam splitter, A/2, A/4 the half- and quarter- 
waveplates, respectively, and A is the iris. 

transition. For lower magnetic fields, the amplitude de- 
pendence is characterized by a single Lorentz resonance 
and the phase is described by an asymmetric curve, both 
centered at twice the Larmor frequency (insets to Fig. [2]) . 
For stronger fields, the single AMOR resonance splits into 
three resonances: the largest central resonance and two 
smaller resonances shifted symmetrically with respect to 
the central one. The splitting is also observed at the cor- 
responding phase dependence. Each of the resonances 
corresponds to different atomic superposition of an indi- 
vidual pair of the ground-state magnetic sublevels. 

As discussed in Sec. UlJ the appearance of the AMOR 
resonance/resonances is associated with generation of the 
ground-state Zeeman coherences. The process is most ef- 
ficient when the modulation frequency matches the fre- 
quency splitting of the magnetic sublevels differing in the 
magnetic quantum number m by 2, Am — 2. In order to 
calculate the splitting of the sublevels, we expand Eq. © 
into the power series of x up to the second order, which 
for the F = 2 state of 87 Rb (/ = 3/2) is equal to 

uj 2 

u m ,m> ~ (m - m!)io L - (to 2 - to' 2 ) L . (20) 

^HF 

For three pairs of Ato = 2 sublevels in the F = 2 state, 
one obtains 

&HF 

w_i,i « 2u L , (21) 
W ,2 « 2oJl +4-r- — . 

&HF 

The first terms in Eqs. (|21l) arise from the linear Zeeman 
effect, while the second ones appear due to the quadratic 
Zeeman effect. For weak magnetic fields, the contribution 
from the nonlinear effect is significantly smaller than the 
ground-state relaxation rate (4w 2 / Ahf < 7). In that 
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Frequency (kHz) 

FIG. 2: Amplitude and phase of the AMOR signal measured 
vs. the modulation frequency. For a magnetic field of about 
640 mG the AMOR signal is split into three resonances due 
to the quadratic Zeeman effect. Each resonance corresponds 
to the different atomic superpositions of the ground-state 
magnetic sublevels. The insets depict unsplit resonance at 
B — 19 mG, where contribution from all superpositions su- 
perimpose. Both signals were measured for a pump power of 
11 fiW tuned to the center of the F — 2 — ► F' — 1 transition 
and a probe power of 3 fiW. 



case, the AMOR resonances associated with individual 
coherences overlap and appear as a single resonance (in- 
sets to Fig. [5]). For stronger magnetic fields, the sublevel- 
splitting frequencies differ sufficiently and the separation 
of the resonances is observed, as seen in Fig. [5J The 
amplitude of the given resonance is determined by the 
amplitude of the corresponding coherence and appropri- 
ate dipole matrix elements. Thus using the NMOR signal 
and also the relations (|17 p -(|19 p determination of the am- 
plitude of the coherence is possible. 

In Fig. [31 the measured splitting of the resonances is 
presented as a function of the magnetic field. For low 
magnetic fields, the resonances are unresolved and no 
splitting is measured [33| . The splitting becomes mea- 
surable for fields corresponding to Larmor frequencies 
above a few tens of kHz and increases quadratically with 
B and uil, which is in good agreement with predictions 
of Eqs. (|2lj). 

In order to verify the model developed in Sec. UH we 
simulate AMOR signals for the first harmonic of the mod- 
ulation frequency, uj m w 920 kHz and a magnetic field of 




Larmor frequency (kHz) 

FIG. 3: (Color online) Average splitting of the AMOR res- 
onances, (uj r — oji)/2, measured vs. the Larmor frequency 
(square). The observed dependences is in very good agree- 
ment with theoretical curve plotted based on Eq. (|21[> (solid 
line) and the data from Ref. |32j]. The data were measured 
with a single sinusoidally modulated light beam of 8-/xW 
power acting as pump and probe simultaneously. 



640 mG (Fig. . The simulations reveal all the features 




_(-)(} I ■ I i i i i I i i i i I i i i i I i i i L 

921.6 922.0 922.4 922.8 

Frequency (kHz) 

FIG. 4: Amplitude and phase of the AMOR signal simulated 
for F = 2 -» F' = 1 transition using Eqs. (JTTJ) . The sig- 
nal reveals all the salient features of real signal (Fig. i.e. 
resonance splitting, similar amplitude relations. 
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observed experimentally (Fig. \2§- For a magnetic field in- 
ducing significant nonlinear Zeeman splitting of the lev- 
els, three AMOR resonances associated with p-2,0, P-1,1, 
and po,2 are observed. The strongest resonance is related 
to the coherence between the | — 1) and |1) sublevels. The 
remaining two resonances are equally split with respect 
to the central one and have equal amplitudes. Also the 
phase of the simulated signal follows the experimentally 
measured dependence crossing zero at the center of the 
largest signal. We attribute the observed deviation of the 
simulated from the measured dependences due to higher 
order processes, such as power broadening and satura- 
tion, which are not included in our model. 

The splitting of the AMOR resonance into three reso- 
nances illustrates a possibility of selective addressing of 
a particular Am = 2 ground-state coherence. Figure [5] 
presents simulations of the density matrix at the first har- 
monic of the modulation frequency calculated at lower 
and higher magnetic fields corresponding to single and 
split resonances (Fig. 2]). The top row is calculated for 
10 mG, whereas the bottom for 640 mG, that corresponds 
to resonance splitting of 120 Hz, (cj r — uii)/2 w 47. The 
simulations are performed for the modulation frequencies 
equal to 2ul ~ 37 (left column), 2ul (middle column), 
and 2wl + 37 (right column). As shown in the top row, 
exact tuning of the modulation frequency to twice the 
Larmor frequency results in generation of all Am = 2 
coherences with the highest and equal efficiency, while 
detuning it away uniformly reduces the amplitudes of all 
of them. In stronger fields the situation is different (bot- 
tom row). While for uj m = 2wl the coherence is 
generated with the highest efficiency, the two other co- 
herences are created significantly less efficient. One can 
selectively increase the amplitude of either of these coher- 
ences by appropriate tuning of the modulation frequency. 
For uj m = 2ujl — 37 (left column), the p-2,0 coherence is 
most efficiently generated, while for uj m = 2u>l + 3j (right 
column), the po,2 coherence has the strongest amplitude. 
This dependence of the amplitude of the specific coher- 
ence on the modulation frequency proves the possibility 
of selective addressing and the control of specific Am = 2 
coherences. 

The AMOR technique is a powerful tool in the anal- 
ysis of a quantum state of the system. In particular, 
scanning the modulation frequency, fitting the data with 
three Lorentz curves, and taking into account strengths of 
the specific transitions, allows one to extract information 
about amplitudes of the density-matrix elements. This 
measurement, however, requires modulation-frequency 
scan within a range strongly exceeding the splitting of 
the resonances. In order not to scan the modulation 
frequency, one may perform free-induction decay mea- 
surements, where information about the coherence am- 
plitude is extracted from time-dependent rotation signal 
(for more details see Ref. [24]|). 

As described above, the theoretical model developed in 
Sec. [II] is valid only for low light powers; for higher light 
intensities such effects as ac Stark shift start to play an 




FIG. 5: (Color online) Calculated absolute values of the 
density-matrix elements (|/9 m , OT '|) of the F — 2 ground un- 
der interaction with linearly-polarized, AM light tuned to the 
F = 2 — > F = 1 transition. The top row corresponds to 
a magnetic field of about 10 mG, when energy splittings of 
all Am = 2 sublevels are equal and the bottom row to a 
much stronger field of 640 mG causing the AMOR-resonance 
splittings of 120 Hz, that is four times the ground-state re- 
laxation rate. The central column shows the density-matrix 
elements for ui m = 2ujl, while the side columns correspond to 
cj m = 2ojl ± 37, respectively. 



important role, e.g. by broadening of the AMOR reso- 
nances. Due to this fact not only the modulation fre- 
quency but also the pump and probe powers determine 
the amplitudes of the Zeeman coherences and hence the 
AMOR signals. In Fig. [6] the amplitude of the AMOR 
signal is presented vs. the pump- and probe-light pow- 
ers. For low pump-light power the AMOR amplitude 
is small (see the upper inset to Fig. 15]), which reflects 
low efficiency of the ground-state coherence generation. 
This efficiency, hence the AMOR-resonance amplitude, 
increases with the pump power and reaches its maximum 
at about 30 fiW. Appearance of the maximum and fur- 
ther decrease of the amplitude results from the higher- 
order effects, such as saturation, hyperfine pumping, and 
repumping/regeneration of the existing Zeeman coher- 
ences. For instance, the hyperfine pumping leads to a 
decrease of a number of atoms in the F = 2 ground state 
transferred to the F = 1 ground state via spontaneous 
emission. Degradation of the AMOR signal is also as- 
sociated with the applied modulation. Efficiency of the 
ground-state coherence generation and hence the number 
of atoms existing in a particular quantum state follows 
the light modulation. At low light intensities, it effects in 
a sinusoidal variation of numbers of atoms evolving with 
specific phases, which results in strong anisotropy of the 
medium. For more intense light, i.e., when saturation 
processes become significant, the efficiency does not re- 
produce the sinusoidal modulation pattern. In particular, 
higher harmonics of the modulation arise in the efficiency 
of coherence generation and number of atoms generated 
during successive pumping phases does not follow the si- 
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FIG. 6: (Color online) AMOR signal as a function of the 
pump- and probe-light power. The data was measured at 
850 mG. The top inset show a cross-section across the plot, i.e. 
amplitude pump-power dependence measured with a probe- 
intensity of 2 /iW. Similarly, the bottom inset shows the data 
taken with varied probe power and fixed pump light intensity 
4 nW. 



nusoidal dependence. It results in weaker anisotropy of 
the medium and a decrease of the AMOR-resonance am- 
plitude for higher pump- light powers (Fig. |6]). 

The dependence of the AMOR signals on the probe- 
light intensity is different. As seen in the lower inset to 
Fig. [6l the amplitude of the AMOR resonance decreases 
with the probe-beam power on the whole accessible range 
of powers. This is caused by the probe light being res- 
onant with the medium (the 80-MHz difference between 
the probe and pump frequencies caused by AOM is neg- 
ligible relative to the Doppler broadening of the tran- 
sition). In such a case, the probe perturbs the atoms; 
absorption of a photon from the probe beam results in 
a new quantum state which, in general, is different from 
the state created initially with the pump-light photon. In 
that way, the probe-light absorption decoheres the sys- 
tem and acts as additional relaxation mechanism reduc- 
ing the AMOR-signal amplitude 

It was shown in Sec. [TT] that efficiency of the ground- 
state coherence generation strongly depends on the pump 
light tuning (dependence on fi). In particular, vari- 
ous pump-power dependences may be observed for the 
light coupling a given ground state with different excited 
states. Such an example is depicted in Fig. which 
shows the ratio of the amplitude of the central resonance 
to the averaged amplitudes of the side resonances vs. 
the pump-light power for the F = 2 —> F' = 1 and 
the F = 2 — > F' = 2 transitions. At low pump pow- 
ers the AMOR signals observed at these two transitions 
are similar in shape and amplitudes, with well-resolved 
triple-component structures [see left-handed side inset to 
Fig.[7£a)]. However, the pump-power dependences of the 
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FIG. 7: (Color online) (a) Ratio of the amplitude of the cen- 
tral of AMOR signal to the averaged amplitude of the side 
components of the signal vs. the light power, (b) ratio of the 
geometric mean of the | — 1) and 1 1} sublevels population to the 
averaged geometric mean of the | —2) and |0) populations and 
|0) and |2) population vs. normalized Rabi frequency. Sig- 
nals were measured at a magnetic field of 790 mG with the 
1.2 ^iW probe power while the calculations were performed 
for a single unmodulated light beam exploiting the rate equa- 
tions. The red points and the red curve corresponds to the 
F = 2 — y F' — 1 excitation while the black squares and black 
curve to the F — 2 — > F' = 2 transition. 



resonances amplitudes measured at the two transitions 
are very distinct. While for the F = 2 F' = 1 transi- 
tion, the ratio increases with the pump-light power [see 
top right-handed inset to Fig. EJa)], and the opposite 
dependence is observed at the other transition. In the 
latter case, the ratio decreases with the power [bottom 
right-handed inset to Fig. EJa)] even below one when the 
side resonances have larger amplitudes than the central 
resonance. The observed dependences reflect the differ- 
ent behavior of the p_n coherence and the p-20 an d P02 
coherences at the two transitions. Among others, the be- 
havior originates from optical pumping and population 
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redistribution between Zeeman sublevels within a given 
hyperfine state. Since for the F = 2 — > F' = 1 transi- 
tion the maximal Clebsch-Gordan coefficients are those 
from the sublevels with maximal m (m = ii* 1 ), these 
states are most efficiently depopulated. Simultaneously, 
the depletion of the other sublevels is smaller, which ef- 
fectively leads to an aligned state with the highest pop- 
ulation in the m — sublevel and the lowest popula- 
tion in the m = ±2 sublevels. The change in the pop- 
ulation distribution is reflected in the amplitudes of the 
Zeeman coherences associated with these states. For the 
F = 2 — » F' = 1 transition and intense light, the lower 
amplitudes have the p-2a and po2 coherences, while the 
higher one has the p~u coherence. The opposite is true 
for the F = 2 — > F' = 2 transition, while the m = state 
is most efficiently depleted. 

In order to qualitatively verify the mechanism de- 
scribed above, the geometric means of the population of 
the magnetic sublevels constituting the coherence were 
calculated vs. the pump intensity. Such a geometric 
mean of the two sublevels' populations sets an upper 
limit on the amplitude of the coherence between the 
sublevels (Schwarz inequality), \p a p\ < ^/PaaPpfi- Fig- 
ure [TJb) shows the ratio of p-\-\p\\ to (y/p-2-2Poo + 
v / pooP22)/2 for the F = 2 state coupled with light to the 
F' = 1,2 states. Populations of the sublevels were cal- 
culated based on the rate equations using CW light and 
neglecting hyperfine optical pumping, saturation, etc. As 
seen in Fig. the simulations qualitatively reproduce the 
observed dependence, i.e. increase of the ratio for the 
F = 2 — > F' = 1 transition and decrease for the other 
transition. 

In alkali atoms there are two hyperfine ground states 
supporting long-living quantum coherences. In particu- 
lar, in 87 Rb there are F = 1 and F — 2 ground states 
separated by about 6.8 GHz that can be selectively ad- 
dressed by appropriate tuning of the pump and probe 
lasers. Figure [5] presents the AMOR signals measured at 
the F = 2 -> F' = 1 and F = 1 -> F' = 1 transitions 
for the same set of experimental parameters. As shown, 
the two signals are significantly different; not only differ- 
ent are the numbers of the AMOR resonances associated 
with the Am = 2 coherences (there is only one Am = 2 
coherence in the F = 1 ground state), but also their 
amplitudes are distinct. The latter difference originates 
from the transition probabilities and less efficient hyper- 
fine pumping at the F = 1 — > F' = 1 transition than at 
the F = 2 —> F' = 1 one. Moreover, the positions of the 
AMOR resonances in a strong magnetic field are differ- 
ent. It results from the difference in the Lande factors 
for the two ground-state 
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where gj is the electron and gi the nuclear (/-factor. 
Based on Eq. (f2"2")l it can be easily shown that although 



FIG. 8: (Color online) AMOR signals recorded for light tuned 
to the F — 2 — > F' — 1 transition (left curve) and the 
F — 1 — i> F 1 — 1 transition (right curve). The difference 
in the position of the AMOR resonances arises from nuclear 
contribution to the Lande factor. Signals were measured at a 
magnetic field of about 780 mG and pump- and probe-light 
powers of 6 pW and 3 pW, respectively. 



the splittings due to the electron spin are opposite and 
hence indistinguishable in the AMOR experiment, the 
nuclear-spin contributions are different, which leads to 
the separation of the resonance Au> m — 2 x 2gip^B /%, 
where p^ denotes the nuclear magneton . For a mag- 
netic field of about 780 mG we predicted a splitting of 
about 4.35 kHz which is consistent with experimentally 
measured value of 4.32 kHz. This difference in frequen- 
cies of AMOR resonances for the transitions is important 
from a point of view of quantum-state engineering since 
it offers an additional way of controlling and modifying 
the quantum state by specific tuning of the frequency of 
an additional RF field j34[. 

Dependences of the ground-state Zeeman-coherence 
lifetimes on the pump and probe powers are presented in 
Fig. [HI The lifetimes r were extracted from the AMOR 
resonance width as r = l/irSu> m , where Suj m is the 
AMOR-resonance half- width at half maximum measured 
vs. the modulation frequency. Figure |H] shows that rais- 
ing the light power of either of the beams leads to broad- 
ening of the AMOR resonance (see insets) and short- 
ening of the ground-state coherence lifetime. In order 
to calculate the coherence lifetime that is not affected 
by the light, we performed a series of measurements of 
the AMOR signals at different pump and probe powers 
and double-extrapolate the resonance width to zero light 
powers. The double-extrapolated lifetime of the coher- 
ences studied in this experiment is equal to 13.2(12) ms, 
which is determined by three relaxation mechanisms: col- 
lisions with the uncoated surfaces, mainly in the cell stem 
containing a rubidium metal droplet, spin-exchange col- 
lisions between rubidium atoms [17| . and temperature- 
dependent dephasing collisions with the cell wall coat- 
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ing [30j. Using a simple mathematical model, the relax- 
ation rate due to collisions with uncoated surfaces was 
estimated at the level of « 2ir x 8.3 s _1 . At the same 
time, the relaxation rate associated with spin-exchange 
collisions, that is calculated based on Ref. [3l|, gives 
2ir x 4.1 s . The other relaxation channel most likely 
are dephasing collisions of atoms with the coating. 



V. SUMMARY AND CONCLUSIONS 

We have analyzed the possibility of generating the 
quantum superpositions of ground-state Zeeman sub- 
levels differing in the magnetic quantum numbers by 2. 
Since in the F > 1 state the sublevels split nonlinearly 
with a magnetic field (nonlinear Zeeman effect), selective 
generation of coherences between specific sublevels by ap- 
propriate tuning of the modulation frequency is possible. 
In particular, it was shown that for the magnetic fields 
such that nonlinear magnetic-sublevel splitting exceeds 
the ground-state relaxation rate, selective addressing of 
the p-20> p-u t and po2 coherences is possible. Efficiency 
of the coherence generation versus different experimental 
parameters, such as modulation frequency, pump- and 
probe-light power, light frequency was analyzed. We 
have shown that in our experimental setup the lifetime 
of the coherences exceeds 10 ms. Such a long coherence 
lifetime opens interesting possibilities for application of 
the coherences in quantum-state engineering and quan- 
tum computation. In this context particularly interesting 
seems to be the ability of modification of atomic quan- 
tum state by application of external field, e.g. magnetic 
and/or RF field. 



FIG. 9: Lifetime of the Zeeman coherence between | — 1) and 
|1) ground-state sublevels vs. pump- (a) and probe- light (b) 
powers. Increasing power of either of the light beams results 
in reduction of the lifetime of a quantum state which manifest 
as broadening of the AMOR resonances (see insets). The 
signals were measured for a magnetic field of about 780 fiG 
with a probe power of 1.2 (a) and pump power of 2.2 fiW. 
The laser was tuned to the center of the Doppler-broadened 
F = 2 -> F' = 1 transition. 
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